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I. IMPORTANCE OF THE GRAVITOMAGNETIC FIELD 



a) A Never-measured Field of Nature 


In electrodynamics, in the frame at rest with an electrically charged sphere 
we have an electric field. If we then rotate the sphere we observe a magnetic field, 

which strength is proportional to the angular velocity. 


Similarly, in Einstein geometrodynamics (Misner et al. 1973, Wheeler 9 ) 

(general relativity), a non-rotating, massive sphere produces the standard 
Schwarzschild field. If we then rotate the sphere we have the occurrence of he 

gravitomagnetic field, whose strength, in the weak field limit, is proportional to the 
sphere angular velocity. In the weak field approximation of the Kerr metric, the 
gravitomagnetic field is (Thorne et al. 1986): 
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where "p = (0,0, - “Y - ) is the gravitomagnetic or Lense-Thirring potential and J is the 
angular momentum of the central body. 

We observe that gravitomagnetism and gravitational waves are the two main 
aspects, still to be measured, of Einstein geometrodynamics analogous to magnetic 
field and electromagnetic waves of electrodynamics. 


b) At the Foundations of Inertia in Einstein General Relativity 


In Einstein geometrodynamics (Misner et al. 1973, Wheeler 1964), to solve the 
initial value problem (York 1979, Choquet-Bruhat and York 1980) as a part of the 
initial conditions, we need to specify on a Cauchy hypersurface L the conforma 


mass-energy current density: j base- Through the initial conditions, with the field 

equations, we then solve for the geometry g ap of the universe and eventually, we 

determine the local inertial frames (where g ap ->n a p) a11 over the spacetime. The axes 

of the local inertial frames (gyroscopes) are therefore influenced, and partially 
determined, by the mass-energy currents in the universe (dragging of inertia 

frames). 


This agrees with some kind of general relativistic formulation of the Mach 
principle, according ro the ideas of Einstein and Wheeler (Wheeler 1988, Car.folini 
and Wheeler 1987). "Nothing would do more to demonstrate the inertia-intiuencing 
effect of mass in motion than to detect and measure the Einstein-Thirring-Lense- 
predicated gravitomagnetism of the earth" (Wheeler 1988). 
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c) A Key Role in Theories of Quasars and Active Galactic Nuclei 

In high-energy astrophysics, some theories of energy storage, power 
generation, jet formation and jet alignment of quasars and active galactic nuclei 
based on the existence of the gravitomagnetic field of a supermassive black 
(Thome et al. 1986). 

In particular, this field may explain the constant direction of the jets over 

millions of light years and therefore over a time of several milhons ; of yea £; 

Through the standard Navier-Stokes equations, in which one includes the 

gravitomagnetic field of the central object, one can show (Bardeen and Petterson 
1975) that the accretion disk tends to be oriented into the equatorial plane of the 
central body - Bardeen-Petterson effect. The jets are then ejected normally to the 
accretion disk, that is, normally to the equatorial plane of the central body The 

angular momentum vector of the central black hole acts, therefore, as a gyroscope 
and this may explain the constant direction of the jets. 

II. LAGEOS III 

LASER RANGED SATELLITES TO DETECT THE GRAVITOMAGNETIC FIELD 
AND SUPPLEMENTARY INCLINATION SATELLITES 
TO AVOID GRAVITY FIELD UNCERTAINTIES 

Many experiments have been proposed to measure the gravitomagnetic field. 

The GPB experiment (Everitt 1974, Lipa et al. 1974) intends to measure the 
Lense-Thirring-Schiff (Lense and Thirring 1918) precession of gyroscopes orbiting 

the earth. 

Polar satellites have been proposed to measure the Lense-Thirnng precession 
of ,hc albLl plane - an onoju/gyroscope (Schiff 1960. Yilmaz 1959) and two 
guided, drag-free, counter-rotating, polar satellites have been suggested to avoid 
orbital inclination errors (Van Patten and Everitt 1976). We recall that a po ar 
satellite: I = 90°, has a null classical nodal precession. 

Here we briefly summarize the new idea (Ciufolini 1985, 1984, 1987, 1988) to 
measure the gravitomagnetic drag of the nodes of two nonpolar, supplementary 
inclination, laser ranged satellites. This idea can be decomposed into two parts: 


( 1 ) 


( 2 ) 


Position measurements of laser ranged satellites, of LAGEOS (Smith and Dunn 
1980 Yoder et al. 1983, Cohen et al. 1985) type, are accurate enough to detect 
the tiny effect due to the gravitomagnetic field: the Lense-Thirring 
precession. 

To cancel out the uncertainties in the enormous perturbations due to the 
nonsphericity of the earth gravity field, we need (Ciufolini 1985, 1984, 1987) 
new satellite, LAGEOS III, with semimajor axis and eccentricity, a and e, equal 
to those of LAGEOS, but with supplementary inclination: I||| = 70 . 

The LAGEOS (Smith and Dunn 1980, Yoder er al. 1985) semimajor axis is 
a = 12270 km, the period P = 3.758 h, the eccentricity e = 0.004 and the inclination 
j _ 109.94°. The period of the node P (Q) = 1046 days. 
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The Lense-Thirring nodal precession is for LAGEOS (Ciufolini 1984). 


. Lense 2GJ ffi 

Q Thirring = , » 2 3/2 = 31 milharcsec/year (2) 

LAGEOS c a (1-e ) 


where J = 5.9 x 10 40 g • cm2/sec = 1.5 x 10 2 cm 2 is the earth angular momentum. 

© 

Unfortunately, the Lense-Thirring precession cannot be extracted from the 
measured value of the LAGEOS nodal precession, because of the uncertainty in the 
theoretical classical precession due to the quadrupole and higher mass moments of 
the earth: 


Classical 

Q = 126°/year 

LAGEOS 

However, a new satellite, of LAGEOS type: LAGEOS-III, with supplementary 
inclination: I, s 70°, would have a classical precession equal in magnitude and 

opposite in sign to that of LAGEOS. By contrast, since independent of the inclination, 
the Lense-Thirring precession (2) would be the same, both in magnitude and sign for 
the two satellites. Therefore, from the sum of the measured nodal precessions, we 

should be able (Ciufolini 1985, 1984, 1987, 1988) to measure the Lense-Thirring effect. 

III. ERROR SOURCES 

ORBITAL INJECTION ERRORS, GRAVITATIONAL AND 
NONGRAVITATION AL PERTURBATIONS AND MEASUREMENT ERRORS 

(1) Errors from gravitational perturbations arise from uncertainties in modeling 
the LAGEOS gravitational nodal perturbations. These errors may be subdivided 
into: 

(a) Errors from orbital injection errors (Ciufolini 1988) due to the 

uncertainties in the knowledge of the static part of the even zonal 

harmonic coefficients, J 2n , of the earth gravity field. These errors are, 

a priori, zero for the two supplementary inclination satellites. However, 
any deviation of the orbital parameters of LAGEOS III from the optimal 
values will introduce uncertainties which must be evaluated. 

(b) Errors from other gravitational perturbations (Ciufolini 1988). Static 

odd zonal harmonic perturbations; static nonzonal harmonic 
perturbations; nonlinear harmonic perturbations; solid and ocean earth 
tides. De Sitter (or geodetic) precession (Bertotti et al. 1987); sun, moon 
and planetary tidal accelerations; nonlinear, n-body, general 
relativistic effects; other very tiny relativistic deviations from geodesic 
motion of LAGEOS. The main error source is, however, due (Ciufolini 

1987, 1988) to the uncertainties in modeling the dynamical part of the 

earth gravity field, that is, to the uncertainties in modeling solid and 
ocean earth tides. 

(2) Errors from nongravitational perturbations (Ciufolini 1987, 1988). Direct solar 
radiation pressure; earth albedo; satellite eclipses; anisotropic thermal 
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radiation; Poynting-Robertson effect; infrared radiation; atmospheric drag, 
solar wind; interplanetary dust; earth magnetic field. The main error sources 

are, however, due to uncertainties in modeling the earth albedo, anisotropic 
thermal radiation, and atmospheric drag (Ciufolini 1987, 1988). 

(3) Errors from measurement uncertainties in the orbital parameters (Ciufolini 
1988). These errors are due to the uncertainties in the measurement of the 
LAGEOS orbital parameters and, in particular, to the errors in the 
measurement of the inclination I and the nodal longitude fl , relative to an 
asympotic inertial frame (Ciufolini 1988). 

IV. PRELIMINARY ERROR ANALYSIS 
A PRELIMINARY 10% MAXIMUM ERROR, 

OVER THE PERIOD OF THE NODE 


A major problem in modeling the LAGEOS orbit is the average secular decrease 

(Cohen et al. 1985, Rubincam 1988) of the LAGEOS semimajor axis of about^ 1 mm per 

day This corresponds to an along-track acceleration of about -3 x 10 cm/sec . 
This acceleration may be explained by three mechanisms (Rubincam 1988, Alonso et 

al. 1988): 

(1) Neutral and charged particle drag (Rubincam 1982, Afonso et al. 1985). 

(2) Thermal drag from the earth infrared radiation (Rubincam, 1987 and 1988) 

' and the thermal lag of the LAGEOS retroreflectors, due to their thermal inertia. 

The re-emission causes an along-track acceleration opposite to the satellite 
motion. 

(3) Thermal drag from the sun radiation plus satellite eclipses by the earth 

(Afonso et al. 1988). The sun radiation is absorbed and then re-emitted by the 
satellite. This phenomenon may be important when the satellite orbits are 
partially in the shadow of the earth. 

The effect of these three perturbations on the node has been investigated in 
relation to the supplementary inclination configuration (Ciufolini 1988, Fannella 

1988). The preliminary result is that their effect should not be larger than 1% 

Q Lense-Thirring due tQ part i c i e drag (Ciufolini 1987, 1988). This is negligible over the 
period of the node for two supplementary inclination satellites due to infrared 
radiation (Ciufolini 1988, Farinella 1988) and not larger than 

2% Cl Lense - Thirrin 8 due to thermal drag plus satellite eclipses (Ciufolini 1988). 


To support these preliminary figures, we observe that to give a secular no a 
precession, a force must: 1) be perpendicular to the orbital plane, and 2) change 

sign in half a period. Even in this worst case, an acceleration with amplitude 

~3 x 1(T 10 cm/sec 2 can, at most, cause a LAGEOS nodal precession of 


5% Cl Lense-Thirring as it j s eas n y calculated (Ciufolini 1988) from the Lagrange 
equation for the node. Ciufolini (1988) has carried out a comprehensive 
preliminary analysis. The result is that, using 2 nonpolar, laser ranged satellites 

with supplementary inclinations, the maximum error, over the period ot the node oi 
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~3 years, should not be larger than -10% of the gravitomagnetic effect to be 
measured. 

A study group, composed of University of Texas at Austin and Italian scientists, 
is, at present, performing a comprehensive computer simulation and covariance 
analysis of the experiment. 


ACKNOWLEDGEMENTS 


D. 

R. 

D. 


I thank C. Alley, E. Amaldi, J. Anderson, P. Bender, N. Cabibbo, M 
Christodoulidis, R. Eanes, F. Everitt, P. Farinella, P. Knocke, L. 
Hellings, H. Mark, R. Matzner, K. Nordtvedt, T. Regge, D. Rubincam, 
Smith, K. Thome, M. Watkins, S. Weinberg, J. A. Wheeler, and C. Will. 


Cerdonio, 
Guerriero, 
I. Shapiro, 


REFERENCES 


Afonso, G„ Barlier, B„ Berger, C„ Mignard F„ and Walch, J. J. 1985, J. Geophys. Res., 90, (Bll), 9381. 

Afonso, G., Barlier, F., Carpino, M„ Farinella, P„ Milani, A., Nobili, A. M. 1988, submitted to A™. Geophys. 
Bardeen, J. M., and Petterson, J. A. 1975, Astrophys. J. Lett., 195, L65. 

Bertotti, B., Ciufolini, I., and Bender, P. L. 1987, Phys. Rev. Lett., 58, 1062. 

Choquet-Bruhat Y.. and York, J. W., Jr. 1980, in: General Relativity and Gravitation, A. Held, ed. (Plenum, 

New York). 

Ciufolini, I. October 1984, Bull Am. Phys. Soc., Cleveland, Ohio 6, 1169, (1985). 

Ciufolini, I. 1986, Phys. Rev. Lett., 56, 278 [submitted October 1984]. 


Ciufolini, I. 1987, Cel. Mech., 40 ,19. 

Ciufolini, I. March 1988, Bull. Am. Phys. Soc.-, Austin, Texas, APS meeting. 

Ciufolini, I. March 1988, A Comprehensive Introduction to the LAGEOS Gravitomagnetic F.xperiment, to e 


Ciufolini I., and Wheeler, J. A. 1987, Geometrodynamics, Inertia Here is Ruled by Mass There (University of 
Texas Press, Austin, in preparation). 

Everitt, C. W. F. 1974, in Experimental Gravitation, B. Bertotti, ed. (Academic Press, New York), p. Si 1. 

Cohen, S. C., King, R. W., Kolenkiewicz, R., Rosen, R. D., and Schutz, B. E., eds. 1985, J. Geophys. Res., 90, 

9215. 


Farinella, P. 1988, to be published. 

Lense, J., and Thirring, H. 1918, Phys. Z., 19, 156. 

Lipa, J. A., Fairbank, W. M„ and Everitt, C. W. F. 1974, in Experimental Gravitation, B. Bertotti, ed. (Arcade 
mic Press, New York), p. 361. 

Misner, C. W., Thorne, K. S., and Wheeler, i. A. 1973, Gravitation (Freeman, San Francisco). 

Rubincam, D. P. 1982, Celest. Mech, 26, 361. 

Rubincam, D. P. 1987, J. Geophys. Res., 92, 1287. 

Rubincam, D. P. 1988, Thermal Drag on LAGEOS, J. Geophys . Res . 

Schiff, L. I. 1960, Phys. Rev. Lett. 4, 215. 

Smith, D. E., and Dunn, P. J. 1980, J. Geophys. Res. Lett., 7, 437. 

Thorne, K. S., Price, R. H., and MacDonald, D. A. 1986, eds., Black Holes, the Membrane Paradigm (Yale 
University Press, New Haven and London), p. 72. 

Van Patten, R. A., and Everitt, C, W. F. 1976, Phys. Rev. Lett. ,36, 629. _ . 

Wheeler, J. A. 1964, Geometrodynamics and the Issue of the Final State, in Relativity, Groups an opo ogy. 

1963 Les Houches Lectures (Gordon and Breach, New York), p. 317. 

Wheeler, J. A. 1988, in IV USSR Seminar on Quantum Geometry, M.A. Markov, ed. (World Scientific, 

Singapore). 

Yilmaz, H. January 1959, Bull Am. Phys. Soc., 4, New York, p. 65. 

Yoder, C. F., Williams, J. G., Dickey, O., Schutz, B. E., Eanes, R. J., and Tapley, B. D. 1983 Nature, 303, 757. 
York, J. W., Jr. 1979, in: Sources of Gravitational Radiation, L. Smarr, ed. (Cambridge Univ. Press, 

Cambridge), p. 83. 


130 



DISCUSSION 


SHAPIRO: LAGEOS-II will of course help to understand the partially unmodeled 

secular variation of the semimajor axis. 


CIUFOLINL Certainly; however this unmodeled variation should not be a problem for 
the gravitomagnetic experiment. In fact, even in the worst case that a comparable 
acceleration would be acting perpendicularly to the orbital plane and changing ig 
in half a period, that is, even in the case of a maximum contribution to the nodal^ rat e 
the total nodal drag due to this acceleration would not be more than 5% of the Lense 
Thirring effect. 


NORDTVEDT: Why do you quote a d J 2/J 2 uncertainty limitation, when you might use 
the LAGEOS orbit, itself, as a measure of J 2 ? 


CIUFOLINI: Unfortunately the nodel precession and the rates of change with time of 

the other orbital parameters are due, not only to the earth quadrupo e moment 
other harmonics too. A solution would be to orbit several high aititude laser-raiiged 
satellites, of LAGEOS type, to measure each even zonal harmonic cocfficie t, h 
proper order, plus one satellite to measure the Lense-Tlurr.ng effect. Another 
solution is to orbit polar satellites (Yilmaz 1959, Eventt and Van Patten 1976), since 

they have a null classical precession. Another solution 18 nf P Qufolini 

inclination satellites to cancel out the classical precession (LAGEOS III, Uutol.m 

1984). 


SHAPIRO: (Comment on question by K. Nordtvedt). The even zonal harmonics (J 2 

J 4 , ...) affect other orbital parameters, as well as W, and these effects for the orbit of 
LAGEOS I, as well as for other satellites, allow estimates to be made of the values of J 2 , 
j 4 It was to the residual errors, or uncertainties, from these estimates that 

Ignazio referred. 
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